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We study the quenched chiral behavior of the pion with mass as low as ~ 180 MeV. The calculation is done on 
a quenched lattice of size 16 3 x 28 and a = 0.2 fm with 80 configurations using overlap fermions and an improved 
gauge action. Using an improved constrained curve fitting technique, we find that the ground state pseudoscalar 
mass versus bare quark mass behavior is well controlled with small statistical errors; this permits a reliable fit of 
the quenched chiral log effects, a determination of the chiral log parameter (S = 0.26(3)), and an estimate of the 
renorrnalized mass of the light quark (m MS {p = 2 GeV) = 3.7(3) MeV). 



1. Simulation Details 

Using a /3 = 2.264 renormalization-group- 
improved Iwasaki 0] gauge action, we study the 
chiral properties of hadrons on a 16 3 x 28 lattice 
with the overlap fermion |^,|| and massive overlap 
operator 



2. Pion Decay Constant and Za 

The renorrnalized pion decay constant is 



f(U) 

J p 



D(mo) = (p + 



m a m a 

-ft) + (p — 7rn^{H) 



where e(H) = H/y/H 2 , H = jsD w , and D w is 
the usual Wilson fermion operator, except with a 
negative mass parameter —p = 1/2k — 4 in which 
k c < k < 0.25; we take k = 0.19 in our calculation 
which corresponds to p — 1.368. 

We use the optimal partial fraction expan- 
sion with a 14th-order Zolotarev approximation 
of the matrix sign function Q]; the sign func- 
tion approximated to better than 3 parts in 10 10 . 
As the conjugate-gradient inverter accommodates 
multi-mass Q we obtain the quark propagator at 
26 masses including 18 masses at or below the 
strange quark mass with less than 10% overhead. 



(ml/2m Q ) (ml/2m ) 
where mo is the unrenormalized quark mass, 
f p U ^ — \/2m- K App is the unrenormalized pseu- 
doscalar decay constant, and App is the am- 
plitude of the two-point local-local correlator, 
G PP (t) — > A P p(e-'"" i + e-'"'( T -')). Since fi R) 
is free of quenched chiral logs, and since we de- 
termine it quite precisely, we use to set the 
scale and obtain a -1 = 0.978(5) GeV. 

As a bonus, we obtain the axial-vector renor- 
malization constant, Za, via the axial Ward 
identity ZaO^A^ — 2Z m moZpP and the rela- 
tions Z^ 1 = Zp = Zs protected by the chiral 
symmetry of the overlap fermion. Thus Za = 

2mp 
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Aj^ which is plotted in Fig. 1. 
A covariant polynomial fit yields 
Z A = 1.85(1) - 0.60(12)m A QCD a 2 

+0.59(4)m 2 a 2 ; A QCD a = 0.25 

so 0(moAQCDa 2 ) and 0(m^a 2 ) terms are small. 
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Figure 1. Za vs. bare mass, mo<x. 

3. Pion Mass and Chiral Logs 

Our plot of (m^a) 2 /tooQ. in Fig. 2, with 18 data 
points at or below the strange quark mass, dra- 
matically reveals the presence of the quenched 
chiral log. We fit m 2 a 2 to the form |9pC[| 
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using a single multi-parameter (weakly) con- 
strained fit (the constraint lifts the degeneracy 
of the three-parameter fit). We obtain simul- 
taneously (the preliminary results) S = 0.26(3) 
and A x a = 1.1(1), which are quite stable against 
changing the window, [m£ un ,m™ a3t ], of masses in 
the fit provided (a) m™ in < 250 MeV; thereafter, 
S monotonically decreases rapidly with increasing 
to™ 11 (e.g. 0.14 for 300 MeV, 0.06 for 400 MeV), 
and (b) m™ ax < 350 MeV; adding a term Bm^a 2 
to the fit formula increases the region of stability 
to TO™ ax > 700 MeV. 

The pseudoscalar decay constant fp has similar 
behavior and fit results to (m 7r a) 2 /mofl. (This is 
automatic for overlap fermions since Z" 1 = Zp = 
Z s implies f P = fi R) V2m 2 j2m .) 
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Figure 2. (m 7r o) 2 /moo vs. bare quark mass, moa. 



Our value of S is significantly higher than other 
recent values []To|— most of which are well below 
the prediction of 6 = 0.183 from the quenched rj' 
loop (using the Witten-Veneziano model of the 
rf mass) in chiral perturbation theory; however, 
we have exact chiral symmetry and, to date, use 
the largest spatial volume and probe the light- 
est pion masses. Indeed, the effect of the chiral 
log makes the data points rise as the pion mass 
drops below to^ w 400-450 MeV. Notice that if 
one restricts the data to pion masses in the range 
of 350-400 MeV to 700 MeV, then the curve can 
be well approximated by a constant within er- 
rors. Equivalently, a plot of m 2 versus too over 
this range would show no deviation from linear- 
ity, as is also seen in 10. The effects of chiral 
logs are seen only at lower pion mass; one is un- 
likely to see the effect at all without data below 
to w w 400 MeV. With adequate spatial volume, 
the extraction of 5 would give a low value with- 
out data below m n w 250 MeV. Our spatial lat- 
tice size is L = 3.2 fm so finite- volume effects are 
expected to be small; our lightest (unrenormal- 
ized) quark mass of w 14 MeV corresponds to a 
pion of about 180 MeV. With a lattice scale of 
a = 0.202(1) fm (set from f„) then L = 3.2 fm 
and (im T ) lightest ~ 2.9. 
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Figure 3. (m^a) 2 vs. pion mass, ra,a. The solid 
curve is the the chiral log fit. The solid lines in- 
dicate the intercepts at the physical values. 

4. Quark Masses 

To renormalize our estimates of the quark 
mass, we follow |j) to determine the renormal- 
ization constants Z^ 1 = Z$ = Zp and thus 
m RGi _ Z m (go)mo(go): The value of the bare 
quark mass mofl that reproduces a given fiducial 
pseudoscalar mass (through x Ie { = (t^to^) 2 , with 
ro/a = 2.885(18) for our lattice) is obtained by 
interpolation of our data for the pion mass. From 
this we obtain Z" 1 = Z$ at three values of quark 
mass (the lowest two quark masses are roughly 
half-strange and strange). The operators and ac- 
tion are 0(a) improved. For Za we find small 
0(77i o AQCDa 2 ) and 0(m 2 a 2 ) dependence as well; 
Za is well approximated (to 1%) by a constant 
up through the strange mass region. Thus we fit 
Z" 1 = Zs to a constant and obtain Z$ = 1.61(7) 
as a first approximation. 

Our lightest pion mass is w 180 MeV so we 
needn't extrapolate far to obtain moa at the phys- 
ical pion mass (Fig. 3). Furthermore, exact chi- 
ral symmetry (zero residual mass) and the form of 
the fit function as determined by quenched ChPT 
provide an end point at zero mass. 



Using m RGI — tuq /Zs and, from the 4-loop cal- 
culation with N f = @, ZIP = Z s /0. 72076, 
we obtain (very preliminary results) for (m u + 



m d )/2: 



,RGI 



5.1(4) MeV and 



MS 



(p = 2 GeV) 



3.7(3) MeV 



Last year's result ]l7j from a 20 4 lattice with 
a Liischer-Weisz gauge action was 4.5(3) MeV. 
(The result quoted in the proceedings was an av- 
erage of this number with an attempt to subtract 
out the effects of the chiral log.) A recent esti- 
mate Q of the world average is 4.5(6)MeV. 
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